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We investigate the reliability of transition strengths computed in the random-phase approximation 
(RPA), comparing with exact results from diagonalization in full Qhuj shell-model spaces. The RPA 
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I. INTRODUCTION 

Electromagnetic and weak transitions provide impor¬ 
tant probes of nuclear wavefunctions. Exact numer¬ 
ical solutions of the many-body Schrodinger equation 
are very difficult, so theorists have devised a number of 
approximations. Transitions offer a serious test of the 
model wavefunctions, not only in comparison to experi¬ 
ment but also between approximations. 

One of the most successful approximations is the in¬ 
teracting shell model (SM), which diagonalizes an effec¬ 
tive Hamiltonian in a restricted space, providing good 
description of the observed low lying states for a wide 
range of mass numbers. The SM wavefunctions suc¬ 
cessfully describe electromagnetic and weak transitions, 
in good agreement (albeit sometimes requiring effective 
couplings) with the experimental data; for this reason, 
in this paper we consider SM results to be exact. Due 
to the computational burden which increases with mass 
number, however, full Ohio SM calculations are presently 
limited to light and medium nuclei {A < 80). 

Another approximation, which works surprisingly well 
considering its simplicity, is mean-field theory, e.g., 
Hartree-Eock (HE) theory. In this approach, the par¬ 
ticles are uncorrelated and move in an average poten¬ 
tial which replaces the two-body interaction; the ground 
state is determined by energy minimization. Starting 
from the HE solution, the Tamm-Dancoff approximation 
(TDA) constructs excited states by assuming them to 
be mixtures of one particle, one hole (Ip-lh) configura¬ 
tions. The HE-I-TDA is the simplest model for ground 
and excited state wavefunctions. There are also limita¬ 
tions of this model: the ground state wavefunction is a 
single Slater determinant which can break symmetries of 
the Hamiltonian and neglects particle-hole correlations 
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arising from the residual interaction P|. 

The random phase approximation (RPA) is a general¬ 
ization of the TDA, which takes into account the residual 
interaction by adding 2p-2h correlations in the ground 
state. (An important extension is the quasiparticle RPA 
or QRPA; because it accounts for nonperturbative pair¬ 
ing effects, it is often considered a superior approach. 
In this paper we restrict ourselves to number-conserving 
RPA.) The excited states, as in the TDA, remain Ip-lh 
mixtures. Is this the case? That is, are the excited states 
mainly mixtures of Ip-lh correlations? The answer is no: 
one knows from SM calculations that, for example, the 
first excited state in has a predominant 4p-4h char¬ 
acter 0, outside the RPA model space. Is therefore no 
surprise that RPA fails in these situations. 

Despite known limitations, the RPA (or QRPA) has 
been widely applied for decades to evaluate spectra and 
transition strengths in nuclei. The reason, aside from 
computational simplicity, is that it is frequently, although 
not always mi, in close agreement with experiment. 
Thus, RPA was the method of choice for description of 
negative parity states in closed shell nuclei iliill 
and in open shell nuclei 1^. RPA and QRPA calcula¬ 
tions (sometimes including continuum states) using phe¬ 
nomenological interactions have been successful in de¬ 
scribing the experimental position of giant resonances 
0, particularly El ma and Ml [13 from electron or 
proton scattering, or Gamow-Teller resonances m . In 
general, however, the description of low-lying transitions 
is poor mui Other studies have used QRPA to com¬ 
pute transitions of interest for astrophysics but did not 
directly compare to experiment [3 • 

What about tests of RPA against exact models? Un¬ 
til recently there was no thorough test of the HE-|-RPA 
against an exact non-trivial model for binding energies. 
We performed such a test for a large number of nuclei 
and found that in general HE-bRPA is a good approxi¬ 
mation to the exact shell-model binding energy [I^, al¬ 
though with some significant failures. Tests of transitions 
have been also mainly against ‘toy’ models [I1II1II3, 



2 


although the QRPA has been tested against exact diago- 
nalization in the full SM space for (3^ or double (3 decays 
[ 23 , m El m , with mixed success. Broader tests of 
RPA transitions strengths against an exact model, e.g., 
SM, has not been done. To address this issue, this pa¬ 
per tests the RPA against full Ofito calculation for several 
nuclei in the sd and pf shells for electromagnetic transi¬ 
tions. Unlike previous tests however, we consider mean- 
field solutions which break the rotational symmetry. 

We find that, in general, the RPA produces a reason¬ 
able approximation to the exact (SM) transitions. The 
most significant failures are for certain low-lying collec¬ 
tive states. We understand these failures through the 
incomplete restoration of broken symmetries and the na¬ 
ture of collective giant resonances. 

The Brown-Bosterli schematic model HEl provides 
insight into giant collective resonances. It assumes the 
model Hamiltonian to be single-particle energies plus a 
separable residual interaction. In both the TDA and the 
RPA, all of the Brown-Bosterli transition strength is to 
a single state, the collective state, which is a model for 
giant resonances. If the residual interaction is repulsive, 
then the collective state will be at high energy. If the 
interaction is attractive, then the collective state will 
be low in energy. In more realistic models, of course, 
the residual interaction includes more complicated two- 
body forces, causing the giant resonance to spread over 
many states. The important lesson of the Brown-Bosterli 
model, however, is that an attractive interaction, such as 
isoscalar quadrupole-quadrupole, leads to large collective 
transitions low in the spectrum, while repulsive interac¬ 
tions, such as cr ■ tj, produce collective transitions lying 
higher in the spectrum. 

Breaking of symmetries can result in low-lying col¬ 
lectivity being subsumed into the ground state. For 
example, the strongly attractive isoscalar quadrupole- 
quadrupole interaction leads to a quadrupole deforma¬ 
tion in the HF state. While the RPA identifies broken 
symmetries, restoration of symmetry is incomplete. This 
point is not fully explicated or understood in the litera¬ 
ture, but our calculations here provide further evidence, 
as discussed in detail in Sec. mill 

The paper is organized as follows. Sec. HD presents 
briefly the interacting SM and the RPA formalisms, with 
emphasis on their application to transition strengths. In 
Sec, mil we compare the SM and RPA transition strengths 
and distribution properties for several nuclides, while in 
Sec.HVI we summarize results. In the Appendix we con¬ 
sider in some detail the effect of breaking of exact sym¬ 
metries on transition strength sum rules. 


II. FORMALISM 
A. General Overview 

In this section we briefly present the treatment of tran¬ 
sitions in both SM and RPA. Before giving any specific 


details, we review some general concepts. 

A primary goal is, given the Hamiltonian H, to calcu¬ 
late solutions to the eigenvalue equation 

H\iy)=E,\i^). (1) 

Suppose we have solved Eq. (^1 , either in the SM or RPA. 
The transition strength from an excited to the ground- 
state is given by the square of the matrix element, 

.S(0^i.) = |HE|0)p, 

while the quantity 

.5fc = ^(F;.-F;o)'=|(HF|0)|2 (2) 

is the (energy weighted) sum rule of order k. In this equa¬ 
tion, V runs over all states, 0 stands for the ground state, 
and F is an arbitrary transition operator. Particularly 
important is Sq, which is the total transition strength 
from the ground state to excited states, and which can 
be rewritten as a ground state expectation value: 

So = J2\{i^\F\0)\^ = {0\F^F\0). (3) 

(In most realistic applications the transition operator is 
a spherical tensor of rank K, Fkm, which has the prop¬ 
erty Then properly the total 

strength in Eq. 0 should be 

So = {Q\Fk-mFkm\Q), (4) 

M 

and similarly for the double commutator in the energy- 
weighted sum rule, etc. To avoid clutter we drop the sum 
over M and it should be assumed to be implicit. ) 

In fact one can write all the sum rules of order k 
as expectation values, most famously the linear energy- 
weighted sum rule: 

= Y.iF, - Fo)mFW)\^ = i(0|[F, [i?,F]]|0). (5) 

We will often characterize our results in terms of 
the centroid, S, and the width, AS”, of the transition 
strengths, defined in terms of Si and S 2 

= <«> 

Both the centroid and the width characterize global prop¬ 
erties of collective excitations. 


B. Shell Model 

In the interacting SM, the number of possible many- 
body configurations is restricted by two means: first, one 
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assumes that only a limited number of nucleons inter¬ 
act (valence particles), the rest forming an inert core, 
and second, the active particles are restricted to a small 
number of single particle states (valence space). Usually, 
the valence space is restricted to a major oscillator shell. 
Diagonalization of an effective Hamiltonian provides the 
low lying states by means of the Lanczos algorithm [^ : 
the corresponding wavefunctions are eigenstates of the 
parity, total angular momentum and isospin. 

The reduction of the available single particle states 
and active particles makes the diagonalization numer¬ 
ically tractable. There is a downside though: in or¬ 
der to take into account the restriction of the Hilbert 
space one has to use modified (effective) operators to 
describe observables or transitions. In many cases a sim¬ 
ple phenomenological use of either enhanced or quenched 
couplings, most famously the former for E2 transitions 
and the latter for Gamow-Teller transitions, yields good 
agreement with experiment. 

In order to calculate transition strengths in SM, we 
have used the Lanczos moment method: the transition 
operator is applied on the initial state wavefunction. 
Then, the newly obtained state is used as the starting 
state, or pivot, for diagonalization by means of the Lanc¬ 
zos algorithm. The size of the pivot vector is the to¬ 
tal transition strength Sq, and the overlap of the pivot 
with the final eigenstate, which it turns out can be read 
off trivially, is the transition amplitude. The interested 
reader is referred to Ref. for details. 


C. RPA 

Several textbooks nm cover the RPA, so we skip 
detailed derivation and just review principal steps. 

The starting point of the RPA is a self-consistent mean- 
field solution, a Slater determinant which can break sym¬ 
metries and which ignores correlations. This determines 
a deformed particle-hole basis, where the occupation 
numbers are zero for particle states and one for hole 
states. 

The RPA ground state is defined as the vacuum for a 
set of quasiboson operators, 

(3,\RPA)=0, 

while the excited states, approximate solutions of Eq. 
o, are given by 

W)=Pl\RPA). 

And because we assume the excited states to be Ip-lh 
correlations, one writes the quasiboson creation operators 

asm 

= E , (7) 

mi 

where 6^^ ^ cl^Ci is the approximate boson mapping 
of the (deformed) fermion operators. We use the con¬ 
ventional notation that m, n are unoccupied (particle) 


states, while i, j are occupied (hole) states, respectively. 
Note that while in Eq. o the first terms describe 
particle-hole correlations on top of the HE state, the 
terms Y^^bmi describe correlations coming from 2p-2h 
configurations in the ground state. In the RPA one as¬ 
sumes that the ground state is still very close to the HE 
solution, so that the hole-particle amplitudes Y are much 
smaller than the particle-hole amplitudes X. Finally, the 
eigenvalue equation 



determines the excitation energies and the particle- 
hole X and hole-particle amplitudes Y. A and B are 
matrices given by 


A, 



c]c„, [H, cl^Ci\ 



(9) 


B 


nj,mt 



[H, clcj],cl^c^ 



( 10 ) 


Thouless showed that if the HE solution corre¬ 
sponds to a minimum in energy surface, the correspond¬ 
ing RPA equation 0 has only real frequencies. In ad¬ 
dition, if the Hartree-Fock state is invariant under some 
particle-hole transformation, such as rotation about an 
axis, this corresponds to a zero frequency RPA mode. 
Thus, the generators of symmetries broken by a mean- 
field solution are eigenvectors of Eq. ® lying at at zero 
excitation energy. This is frequently interpreted as “ap¬ 
proximate restoration of broken symmetries” Q; in fact, 
it is more accurate to say that the RPA respects sym¬ 
metries by separating out exactly spurious motion. In 
Sec. nni we present evidence that significant part of the 
quadrupole response for even-even nuclei is still retained 
in the ground state, which suggests that broken symme¬ 
tries are only partially restored by the RPA. 

The RPA provides a model for excited states, and 
to calculate the transition probability from any non- 
spurious state to the ground state one needs the tran¬ 
sition matrix element (;/|F|RPA). In the RPA, the latter 
can be written in terms of particle-hole amplitudes X'^ 
and Y'^, namely, if one has a one-body transition opera¬ 
tor (and a spherical tensor of rank K), Fkm, which can 
be written as 


Fkm = '^{fZ4nC^ + (-l)^/„^^^4cm), (11) 

mi 


then Pi 

(j.|Ekm|RPA) = /“ . A" + (-1)“/-^ . Y-' (12) 

where f ■ X = fmiXmi etc. 

With the transition matrix element it is possible 
to calculate in the RPA any moment of the distribution 
strength, and therefore the total strength, the centroid 
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and the width. Sec. IIIII compares the SM and RPA pre¬ 
dictions for these quantities, as well as individual transi¬ 
tion strengths in several nuclei. Before proceeding with 
our numerical results, there are two points to discuss. 

First, the energy-wei ghte d sum rule ^i. The RPA has 
the famous property [J, : 

^f7.1(z.|i^|RPA)p = i(HF|[F,[if,F]]|HF}. (13) 

V 

In the Appendix we revisit the derivation of Eq. dia and 
find that it can be violated */ an exact symmetry such as 
rotational invariance is broken. In Sec. cudi we confirm 
the violation numerically, and find the worse case to be 
where the bulk of the transition strength lies very low 
in energy, such as isoscalar E2. This appears to be a 
new result. Because applications of RPA have usually 
assumed spherical symmetry they are not invalidated, 
but ambitious RPA calculations that allowed for broken 
symmetries in the mean-field [31)) should be approached 
with caution. Of particular interest, which we have not 
yet explored, is the QRPA which breaks particle number 
conservation. 

A second point we would like to discuss involves 
ground-state to ground-state transitions. (These are 
nothing more than ground-state expectation values of 
non-scalar operators. We discuss ground state expec¬ 
tation values of scalar operators in |3ll| : in principle one 
can extend such calculations to non-scalar operators, but 
we have not yet done so. Marshalek and Weneser 
discuss expectation values for electric quadrupole and 
magnetic moments, but their approach is not very trans¬ 
parent for general implementation.) Most discussions of 
the RPA do not give a well-defined procedure to calcu¬ 
late such transitions, in part because they vanish when 
the Hartree-Fock state has spherical symmetry. While 
such g.s. to g.s. transitions are forbidden for a spherical 
(hence J = 0) state, they are in general not forbidden 
for nonspherical HF states. If RPA does not fully restore 
broken symmetries, a significant contribution to the to¬ 
tal strength could be absorbed into otherwise forbidden 
g.s. to g.s. transitions. We investigate in detail this point 
in Sec. irrm and indeed we find that significant strength 
to excited states can be missing for even-even nuclides. 


III. RPA VS. EXACT SHELL-MODEL 
STRENGTHS 

In order to test the RPA reliability for computing tran¬ 
sition strengths, we calculate both the mean-field and 
exact solutions in the same model space, using the same 
Hamiltonian. 


A. Model space, interaction and transition 
operators 

We work in full Ohw shell model spaces, restricting 
the single particle states to one major shell. Most of 
our examples were computed in the sd shell, limiting 
the nucleons outside an inert core to the single 
particle orbits lsi/2-0fi3/2-0(i5/2. Additionally, we con¬ 
sidered two nuclei, “^^Ti and in the pf shell, i.e., 
lPi /2 ~ lP 3/2 ~ 0 / 5/2 “ 0 / 7/2 single particle states outside 
“^^Ca core. For the interaction, we used the Wildenthal 
“USD” in the sd shell |3^ and the monopole-modified 
Kuo-Brown “KB3” in the pf shell We emphasize 
that due to our restriction to a single major shell and lim¬ 
itation to mixing angular degrees of freedom, the mean- 
field solution can break only the rotational symmetry. 

For testing purposes, we have considered Fjt = 
err'^Yj, with J = 2, that is isoscalar (T = 0) and 
isovector (T = 1) electric quadrupole (E2). While good 
agreement with experimental transitions strengths re¬ 
quires nontrivial effective proton and neutron charges, 
the main contribution of effective charges is a rescaling 
of the strengths; therefore for simplicity we took the bare 
charges, = 1 and e„ = 0. This might appear to sug¬ 
gest that we only considered the proton response, but 
we computed the isoscalar and isovector responses sep¬ 
arately, that is, incoherently (only if the isoscalar and 
isovector responses were summed coherently would one 
have pure proton response). 

In addition, we tested transition distributions for spin 
flip (SF) and Gamow-Teller (GT) which are the isoscalar 
and isovector components of the spin operator a. To 
avoid confusion, note that the actual GT operator we 
used is (JTz so that = Z — N was conserved. 

A large fraction application of RPA calculations are to 
El transitions. Because our shell model valence space 
does not include single-particle states of opposite parity, 
we could not investigate El transitions here. 


B. Results for isovector quadrupole, SF, and GT 
transition operators 

In this section we show results for isovector E2, SF 
(spin-flip) and GT (Gamow-Teller) transition operators. 
The main common feature is that their collective tran¬ 
sitions lie relatively high in energy. We find that for 
such transitions the RPA is in reasonably good agree¬ 
ment with the SM results, especially for the total transi¬ 
tion strength. 

Figures HE compare the RPA and SM transition 
strengths; we choose for exemplification ^°Ne (even- 
even), ^^Ne (even-odd) and ^^Na (odd-odd), but the gen¬ 
eral trend is the same for all the nuclides investigated. 
The excitation spectra are discrete, but to guide the eye 
we folded in a Gaussian of width 0.7 MeV. In addition, 
tables nnni summarize the results in both SM and RPA 
for several nuclei; we present only the total strengths, the 
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TABLE I: Total strength So, centroid S, and width AS 
for isovector E2 transition operator. The nuclei have been 
grouped into even-even, odd-odd and odd-A. 


So_ S (MeV) AS (MeV) 


Nucleus 

SM 

RPA 

SM 

RPA 

SM 

RPA 

2°Ne 

0.98 

1.15 

14.53 

11.92 

3.47 

2.44 

22 Ne 

2.37 

1.86 

8.15 

7.70 

5.85 

4.25 

'"Mg 

1.88 

1.96 

14.40 

11.86 

4.09 

2.46 

28Si 

2.28 

1.96 

14.35 

13.41 

4.29 

1.93 

3«Ar 

1.38 

1.34 

12.49 

11.01 

4.24 

3.56 

44 Ti 

2.15 

1.90 

8.23 

6.68 

2.86 

1.98 

22 Na 

1.60 

1.69 

11.67 

10.44 

4.37 

2.61 

24Na 

2.07 

2.09 

9.82 

8.21 

6.26 

4.14 

46y 

2.32 

3.00 

7.96 

6.62 

4.17 

1.97 

21 Ne 

1.39 

1.43 

10.52 

9.41 

5.66 

4.27 

25Mg 

2.28 

2.20 

11.48 

9.47 

6.21 

4.41 

2®Si 

2.52 

2.20 

11.61 

10.21 

5.68 

4.25 


centroids and the widths of the distributions. 

The figures show that the RPA calculations follow the 
general features of the SM transition strength distribu¬ 
tions. Note however that by comparison to SM, the RPA 
distributions have smaller widths (see tables HHIIll . This 
is not surprising, as higher-order particle-hole correla¬ 
tions are expected to further fragment the distribution. 
The RPA centroids are generally shifted to lower energies 
than the SM. Although the centroids are related to the 
energy-weighted sum rule Si , we remind the reader that 
we do not violate Eq. because the HE state is only 
an approximation to the ground state. Furthermore, the 
shift in the centroid does not appear correlated with the 
correctness of the RPA estimations of the ground state 
energy or other observables [sj. One might expect 
that the correct inclusion of the pairing interaction by 
means of HFB-I-QRPA would improve the results. This 
is reasonable and worth trying, but see discussion and 
caveats regarding pairing and QRPA in 

For computational simplicity, we restrict ourselves to 
real wavefunctions; this has no effect for even-even nu¬ 
clei. But because the rotations about x or z axis are 
complex, for odd-odd or odd-A nuclei the RPA does not 
identify all the corresponding generators as exactly zero- 
frequency modes. Instead, we obtain a ‘soft’ mode at very 
low excitation energy. Transition strengths to the soft 
mode are in fact ground-state-to-ground-state strength 
normally not computed in RPA. 

To summarize the results in this section, we have com¬ 
pared the SM and RPA distribution strengths for isovec¬ 
tor E2, SF and GT transition operators. We found in 
general good agreement for the total strength in several 
nuclei. While less satisfactory, the centroids and widths 
of the distributions are still close. As a general feature 
however, the RPA distributions are smaller in width and 
lower in energy than the SM results. 



FIG. 1: Isovector E2, SF and GT transition strengths for 
^°Ne. Both the exact SM (solid curve) and RPA (dashed 
curve) distributions have been smoothed with a Gaussian of 
width 0.7 MeV to facilitate comparison. 



FIG. 2: Same as in Fig. Qfor ’^^Ne. 
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Ex. Energy (MeV) 

FIG. 3: Same as in Fig. 0for ^^Na. 

TABLE II: Same as in Table |I] for SF transition operator. 


So_ S (MeV) AS (MeV) 


Nucleus 

SM 

RPA 

SM 

RPA 

SM 

RPA 

2°Ne 

1.05 

1.23 

17.10 

12.26 

4.38 

1.95 

Ne 

3.53 

4.44 

11.40 

8.82 

4.32 

2.37 

'"Mg 

4.15 

4.78 

13.22 

10.17 

4.48 

1.97 

CO 

5.82 

5.20 

12.75 

11.62 

4.34 

1.89 

36 Ar 

2.68 

2.70 

14.53 

11.17 

3.69 

2.88 

44 Ti 

2.56 

3.32 

9.98 

7.86 

2.56 

1.55 

22 Na 

8.57 

5.78 

5.01 

6.67 

5.22 

3.49 

24Na 

10.06 

7.66 

5.83 

7.05 

5.48 

3.19 

46y 

5.44 

7.68 

8.76 

6.40 

2.51 

2.34 

24Ne 

4.02 

3.54 

7.50 

7.62 

5.92 

3.95 

25Mg 

6.94 

6.33 

9.19 

8.47 

5.73 

3.36 

29Si 

8.42 

8.47 

9.38 

7.86 

5.07 

4.45 


C. Results for isoscalar quadrupole response 

This section presents comparison between the SM and 
RPA distribution strength for the isoscalar quadrupole 
transition operator. The main difference with respect to 
the other transitions investigated in this paper is that 
the collective strength lies very low in energy, for re¬ 
alistic Hamiltonians have a strong attractive isoscalar 
quadrupole-quadrupole component. 

We considered again for comparison the same nuclides 


TABLE III: Same as in Table ||]for GT transition operator. 


So_ S (MeV) AS (MeV) 


Nucleus 

SM 

RPA 

SM 

RPA 

SM 

RPA 

20Ne 

1.05 

1.33 

16.32 

12.53 

4.35 

2.42 

^^Ne 

3.87 

4.85 

12.00 

9.37 

4.48 

3.16 

'"Mg 

4.26 

4.85 

14.46 

11.74 

4.24 

2.42 

3®Si 

6.65 

5.70 

15.19 

13.77 

3.59 

1.88 

36 Ar 

2.74 

2.79 

14.85 

12.09 

3.45 

2.99 

44 Ti 

3.03 

3.74 

10.12 

8.42 

2.86 

2.43 

32Na 

5.51 

5.47 

9.96 

9.28 

4.35 

3.18 

’24Na 

7.43 

7.71 

10.32 

9.29 

4.87 

3.48 

46v 

10.60 

7.85 

4.93 

8.15 

4.37 

2.28 

34Ne 

4.25 

3.55 

7.87 

8.67 

5.97 

3.98 

26Mg 

7.12 

6.76 

11.02 

10.00 

6.05 

4.21 

36Si 

9.42 

8.63 

12.28 

10.39 

5.41 

4.99 


> 



FIG. 4: Isoscalar E2 transition strengths for ^°Ne. The SM 
(solid curve) and RPA (dashed curve) distributions have been 
smoothed with a Gaussian of width 0.7 MeV. The large col¬ 
lective peak at low but nonzero excitation energy for the SM 
is absent in the RPA; see text for discussion. 


investigated previously, and we plot the SM and RPA 
distributions in figures HE Characteristics of the dis¬ 
tributions for several other nuclei are given in table HU 
In contrast with the results in Sec. IIII Bl we End a large 
discrepancy between the total strengths in RPA and SM, 
especially for even-even nuclei. 

Figure^lshows that, if one ignores the low energy tran¬ 
sitions, one obtains again a reasonable agreement be¬ 
tween the SM and RPA distributions. Similar features 
encountered for other transitions appear, that is a lower 
energy centroid and smaller width of the RPA distribu¬ 
tion with respect to SM. 

As for the relative good agreement for odd-odd and 
odd-A nuclei, we have to point out that most of the RPA 
strength is concentrated in the lowest energy state which, 
as already noted, appears just as an artifact of our ap¬ 
proach (restriction to real numbers). A full treatment 
of rotations by inclusion of complex numbers would shift 
these ‘soft-mode’ states to zero modes, that is, degenerate 
with respect to the ground-state, and we would expect 
the odd-odd and odd-A cases to then resemble the even- 
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FIG. 5: Same as Fig. glfor ^^Ne. 


Isoscalar E2 
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FIG. 6: Same as Fig. glfor ^^Na. 


TABLE IV: Same 
operator. 

as in Table 0 for isoscalar E2 transition 

Nucleus 

50 

S (MeV) 

A5 (MeV) 

SM 

RPA 

SM 

RPA 

SM 

RPA 

2°Ne 

7.86 

0.19 

2.12 

9.81 

1.92 

2.30 

^^Ne 

9.36 

0.89 

2.01 

5.52 

2.19 

2.79 

'"Mg 

12.57 

0.51 

2.13 

7.99 

2.09 

2.75 

"«Si 

12.04 

0.56 

2.51 

9.88 

2.33 

2.33 

36 Ar 

7.17 

0.23 

2.42 

9.57 

1.91 

2.74 

44 Ti 

10.87 

1.50 

1.73 

3.99 

1.73 

1.70 

32 Na 

9.53 

7.49 

1.47 

1.27 

2.63 

1.82 

34Na 

8.81 

6.33 

2.10 

1.81 

2.85 

1.88 

46v 

15.21 

15.20 

1.62 

0.87 

1.94 

1.63 

34Ne 

8.74 

13.27 

1.53 

0.64 

2.82 

1.35 

33Mg 

10.71 

12.49 

2.25 

1.08 

2.66 

1.62 

3®Si 

9.70 

1.38 

2.72 

4.66 

2.62 

4.25 


even cases: missing the low-energy collective strength. 
(Note that qualitatively the results for ^®Si, for which 
we obtain the correct number of zero RPA modes, are 
similar to the even-even nuclei.) Conversely, we can turn 
around these results into a hypothesis: that the missing 
low-lying collective strength in even-even nuclides are due 
to incomplete symmetry restoration, and that the miss¬ 
ing strength resides in the RPA ground state. The fact 
that the missing strength shows up in soft modes that 
arise as artifacts of our computational methods bolsters 
this hypothesis. For the interested reader, more details 
can be found in the next section. 


D. “Missing strength” and broken symmetries 

In this section we provide further evidence supporting 
our hypothesis that the low-lying collective strength is 
missing due to incomplete symmetry restoration in the 
RPA, and a significant fraction of the RPA strength gets 
absorbed in a ground state to ground state transition. 

Our first test of incomplete symmetry restoration is the 
comparison of the transition strength for spherical and 
deformed HF solutions. While the proton-neutron inter¬ 



FIG. 7: ^®Si: Isoscalar E2 for deformed and spherical HF 
state. 


action induces deformation in the HF Slater determinant 
for ^®Si, it is possible nevertheless to force a transition to 
a spherical HF state (both protons and neutrons filling 
the d^ji orbits) by increasing the gap between the d^ji 
single particle energy and the other single-particle states. 
The above-reported values for ^®Si use the USD value 
of e((i 5 / 2 ) = —3.94 MeV, which yields a deformed HF 
state. In addition we computed ^®Si at €((( 5 / 2 ) = —5.64 
MeV and —5.74 MeV. At —5.64 MeV the HF state is 
still deformed while at —5.74 MeV the HF state is spher¬ 
ical. (Actually, for these values of £((( 5 / 2 ) there exist both 
spherical and deformed locally stable HF solutions, but 
at —5.64 MeV the deformed state has a slightly lower 
HF energy while at —5.74 MeV the spherical state has 
the lowest HF energy. Thus this is a first order ‘phase 
transition’ as described in section 4 of . The so-called 
collapse or breakdown of RPA, readily seen in toy models 
such as the Lipkin model 0, only occurs when one has 
a second order ‘phase transition,’ when one has only one 
stable HF solution. In other words, our RPA calculations 
do not collapse at the transition point.) 

Fignre[7|shows small difference in the SM strength dis- 




















tribution in contrast with a dramatic change for RPA. 
The difference between the (^ 5/2 single-particle energies 
in the two cases is small and one can follow a smooth 
change for all observables in the SM; we have therefore no 
reason to suspect any fundamental difference in the struc¬ 
ture of the states. Note however that, when the HF state 
is spherical, the low lying states are correctly described 
in the RPA, the reason why the RPA was successful in 
describing low lying collectivity in closed shell nuclei; but 
note also that the high-lying part of the strength is not 
correctly described. In contrast, when the HF state is de¬ 
formed, the RPA strength distribution changes dramat¬ 
ically, even through the SM strength distribution does 
not: the low-lying strength vanishes, but the high-lying 
strength is approximately correct. 

As a second test, we compare the total strength and 
the energy-weighted sum rule Si computed in different 
ways. 

Table El presents the total strength Sq for a transi¬ 
tion operator F, where F is either the spin-flip operator 
or isoscalar E2 operator. The columns labled ‘SM’ are 
the exact shell model results, for which Sq = ( 0 |F^| 0 ) = 
l(^l-^| 0 )P- Of course, for the shell model both meth¬ 
ods yield the same result. The columns RPA-X and RPA- 
E correspond to equivalent methods for the RPA. RPA-X 
is the expectation value {RPAlF^^F\RPA) as laid out in 
1^. where we showed that the RPA expectation value 
was often a reasonable approximation to the shell model 
result, though not always as seen for some of the spin- 
flip cases. RPA-E is the sum where 

the sum is only to excited states as g.s.-to-g.s. transi¬ 
tions are difficult to define in RPA. 

The horizontal lines in TablelVisegregate even-even nu¬ 
clei with deformed HF states, even-even spherical nuclei, 
odd-odd nuclei, and odd-A nuclei. Keep in mind that 
for the latter two groups we do not get all the true zero 
modes (because we restricted the Slater determinant to 
real single-particle wavefunctions), but at least one zero 
mode is replaced by a soft mode, except for the case of 
^®Si which does have all the expected zero modes. 

What do we learn from TableEf We draw the reader’s 
attention to the isoscalar E2 strength in deformed even- 
even-nuclides, and in ^®Si, all of which have the ex¬ 
pected number of exact zero modes. Here the summed 
RPA strength (RPA-E) is dramatically and consistently 
smaller than either the exact SM result, or the expecta¬ 
tion value RPA-X. By way of contrast, the nuclides with 
spherical HF states and thus no zero modes, or those who 
have soft modes rather than zero modes, have summed 
RPA strength in reasonable accord with the SM total 
strength. Furthermore the RPA expectation value oiQ-Q 
also agrees with the SM total strength, which suggests to 
us that some of the missing RPA strength is in g.s. to 
g.s. transitions. This line of reasoning is weakened by 
the poor reliability of the RPA expectation value, as dis¬ 
cussed in and as seen in the spin-flip values, which 
take on unphysical negative values for spherical nuclei. 

Therefore, to further dissect this issue, in Table 


TABLE V: Comparison total strength So as computed in the 
shell model (SM), by taking the expectation value of an ob¬ 
servable (RPA-X), and by summing the RPA strengths di¬ 
rectly (RPA-E); see text for more details. The horizontal 
lines separate even-even deformed nuclei, even-even spheri¬ 
cal nuclei, odd-odd and odd-A. Notation on ^®Si: ^ indicates 
£d 5/2 = —5.74 MeV and a spherical HF state. 




SF 



Isoscalar E2 

Nucleus 

SM 

RPA-X 

RPA-E 

SM 

RPA-X 

RPA-E 

2°Ne 

1.05 

1.34 

1.23 

7.86 

8.17 

0.19 

’^^Ne 

3.52 

1.94 

4.44 

9.36 

9.98 

0.89 

""Mg 

4.15 

5.53 

4.78 

12.57 

12.52 

0.51 

36 Ar 

2.68 

2.90 

2.70 

7.17 

7.66 

0.23 


5.82 

5.13 

5.20 

12.04 

13.77 

0.56 

28git 

9.07 

-5.22 

12.53 

7.95 

7.50 

6.86 

22o 

5.04 

-0.60 

6.31 

2.71 

2.75 

2.36 

24o 

5.17 

-1.07 

6.24 

1.88 

1.92 

1.57 


8.57 

7.64 

5.78 

9.53 

9.92 

7.49 

^■‘Na 

10.06 

6.97 

7.66 

8.81 

9.61 

6.33 

2iNe 

4.02 

1.77 

3.54 

8.74 

8.69 

13.27 

^^Mg 

6.94 

1.52 

6.33 

10.71 

12.34 

12.49 


8.42 

5.77 

8.47 

9.70 

11.67 

1.38 


ED we consider the energy-weighted sum rule Si. 
The SM value is “ Eo)\{’y\F\0)\^. The RPA 

value is f2,y|(iz|E|RPA)p, while the HF value is 
^(HFj)!^, [i7, EJIHF) (for technical reasons, discussed in 
the Appendix, we can only compute this for even-even 
nuclei). 

For nuclei with spherical HF states, that is, no zero 
modes, the RPA Si and the HF value are identical; this 
is the usual theorem regarding the energy weighted-sum 
rule. For nuclei with deformed HF states, and thus 
with zero modes, the RPA and HF values differ, a small 
amount for isovector E2 and dramatically for isoscalar 
E2. (We also computed but do not show values for the 
spin-flip operator; the results are similar to isovector E2, 
although with an even smaller contribution from the zero 
modes.) Overall these results are consistent with our 
hypothesis that low-lying strength is subsumed into the 
RPA ground state due to incompletely restoration of a 
broken symmetry; the difference is larger for isoscalar E2 
because of the large strength at low energy. 


IV. CONCLUSIONS 

The purpose of this paper was to investigate the reli¬ 
ability of the RPA for calculating transition strengths in 
nuclei. To accomplish this we have computed the RPA 
and shell model strength distributions in the same Ohio 
shell model space. 

The comparison between RPA and SM showed two dif¬ 
ferent results, depending upon the nature of transitions. 
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TABLE VI: Comparison of energy-weighted sum rule as 
computed in the shell model (SM), taking the weighted sum of 
RPA strengths (RPA), and taking the HF expectation value 
of a double-commutator (HF); see text for more details. No¬ 
tation on ^®Si, same as for Table El 


Isovector E2 Isoscalar E2 


Nucleus 

SM 

RPA 

HF 

SM 

RPA 

HF 

2°Ne 

14.27 

13.74 

13.74 

16.63 

1.82 

7.43 

Ne 

19.28 

14.34 

14.40 

18.84 

4.92 

10.51 

2^Ne 

21.89 

14.89 

15.06 

20.87 

8.42 

11.99 

'"Mg 

27.09 

23.28 

23.28 

26.71 

4.08 

14.87 

36 Ar 

17.24 

14.72 

14.72 

17.32 

2.21 

8.64 

26Si 

32.66 

26.22 

26.22 

30.22 

5.58 

17.67 

28sit 

40.13 

35.76 

35.76 

34.31 

28.26 

28.26 

220 

11.46 

8.56 

8.56 

11.46 

8.56 

8.56 

220 

10.42 

7.99 

7.99 

10.42 

7.99 

7.99 


Thus, we found that when the strong collectivity lies at 
high energies, such as isovector E2, SF and GT tran¬ 
sitions, the RPA and SM are in reasonable agreement. 
When the transitions lie at low energies however, the 
agreement is poor. We presented evidence that the prob¬ 
lem arises from an incomplete restoration of the sym¬ 
metries broken by the mean-field; for low-lying transi¬ 
tions we propose that significant part of the transition 
strength is subsumed into the RPA ground state. Future 
work should directly investigate ground-state to ground- 
state transitions in the RPA. (These are also needed for 
ground state moments, such as magnetic dipole or elec¬ 
tric quadrupole, of odd-A nuclides.) Finally, we also have 
found, both analytically and numerically, that the stan¬ 
dard lore regarding the RPA energy-weighted sum rule 
Fq. (II dll does not hold if an exact symmetry is broken, 
particularly if the centroid of the transition strength is 
very low in energy. 

This paper also marks a final stage within a larger 
project to test the reliability of the HF-I-RPA for a 
global microscopic theory of nuclear properties min 
We conclude that HF-I-RPA is a good starting point 
for such a task, but because of occasional failures fu¬ 
ture work should investigate Hartree-Fock-Bogoliubov 
-l-QRPA and extensions such as renormalized RPA, self- 
consistent R PA, etc., (see Ref. P| as well as the bibliogra¬ 
phies of 1^) and the second RPA, which has been 
shown to differ significantly from the standard RPA in 
its description of F2 giant resonances of Our 

work suggests an important and specific test of any such 
“improvement” to RPA: the description of low-lying col¬ 
lective strength, such as isoscalar F2, which is sensitive 
to restoration of the rotational symmetry broken by a 
deformed mean-field state. 


APPENDIX: THE ENERGY-WEIGHTED SUM 
RULE, REVISITED 

In Table EH we saw a discrepancy between two ways 
to compute the energy-weighted sum rule, 

= i^(-l)'^(HF|[FK.-M,[i?,FKM]]|HF) (A.l) 

M 

and 

n,\{i^\FKM\RPA)\\ (A.2) 

M 

In textbooks 0 one finds “proof” that Si^ = 
which is often stated that RPA preserves the energy- 
weighed sum rule. In this Appendix we revisit the 
proof, with special attention to zero modes that arise 
from broken exact symmetries, and we find that instead 
a term that arises from zero modes. 

Suppose we have a broken symmetry, such as rotational 
invariance. The Hartree-Fock state is deformed and has 
a particular orientation, but the Hartree-Fock energy is 
independent of the orientation. This shows up in the 
RPA matrix equation 0 as a zero-frequency mode. For 
H > 0 one has the normalization = 1 but this 

normalization is impossible for H = 0. Instead one intro¬ 
duces collective coordinates Qi, and conjugate momenta 
which satisfy 

aA - BP* = iMMlQ., 

AQ,-BQ: = -^P^. (A.3) 

Here M^, is a constant, interpretable as mass or moment 
of inertia fixed by the normalization of P,Q'. 

Q1-P.-Qx-P: =iSx.. (A.4) 

Note that if A and B are real, then X, Y are real, but of 
necessity P and Q are complex (one is real and the other 
imaginary). With these zero-mode frequencies one must 
supplement the quasiboson operators /3, in Eq. o 
with the generalized coordinate and momentum opera¬ 
tors Q,V. 

Because of the expansion CD one can use the defini¬ 
tions ®, 113, and able to use A and B to write 
as 

^ ^ (A.5) 

M mi,nj 

The matrices A and B can be written in terms of par¬ 
ticle hole amplitudes X and Y, and the canonical mo¬ 
mentum operators associated with broken symmetries 

A = XnX^+ Y*PLY'^ + PM-^P\ (A.6) 
-B = XVlY'^+Y*VlX^ - PM-^P'^. (A.7) 
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Substitution and some algebra yields 


gHF _ gRPA 




-Pf. 


(A.8) 


(-r/ 


M f-M 




Although we do not show it, one can write the right-hand 
side as a double-commutator of boson operators. 

As further motivation, one can start from Eq. 03 and 
write the contribution from a single frequency to the RPA 
energy-weighted sum rule as 


that goes into g.s.-to-g.s. transitions, due to incomplete 
symmetry restoration. Undoubtedly more work remains, 
but we hope our results act to inspire further careful in¬ 
vestigation. For example, for some transition operators 
there is no or very small contribution from the zero modes 
even for nuclides with deformed HE states; this seems to 
be associated with transitions with high-lying giant reso¬ 
nances, again consistent with our interpretation of incom¬ 
plete restoration of symmetries and low-lying strength 
being subsumed into the RPA ground state. 

Finally, we wish to discuss the general rule 


M 


Even if g. is not a zero mode, one is free to transform to 
collective coordinates and momenta Q, 


= 


yn — _ 






(A. 10) 






P* (A.ll) 


Inserting into (I A. Oil and letting ^ 0, there is a fi¬ 
nite remainder exactly equal to the rightmost term of 
(Toil . It is of course surprising to find contribution to the 
energy-weighted sum rule from ‘zero excitation energy’. 
But we remind the reader that the RPA zero modes are 
continuous states arising from broken symmetries in the 
mean-field solution, while the discrete states are quan¬ 
tized harmonic oscillations about the HE Slater determi¬ 
nant. Part of the challenge of RPA is that continuous 
modes clearly cannot be treated on the same footing as 
the discrete states. 

We find numerically that the discrepancy in Table ED 
is given exactly by the last term in Eq. lA.QII . Our inter¬ 
pretation of Table EH and Eq. 1A.9II is missing strength 


- E^)MF\r)\^ = i(M|[F, [H,F]]\^l). (A.12) 


Eq. IXn is true for any true eigenstate |/r), and so holds 
for full shell-model calculations. But the HE state is not 
an eigenstate, and so one cannot use the Lanczos moment 
method described in Section EH Instead, we take the 
HE state and project onto a vector in the basis of shell- 
model states (because shell-model basis states have good 
Jz, the projection can only be done easily for even-even 
nuclides); both FI and Fkm are matrices in the restricted 
model space and we compute directly {HF — Fi?)|HF) 
and dot that vector onto E|HF). One must sum over 
all shell-model states with intermediate values of Jz, a 
tedious but necessary task for computing . 
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